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Wc report on a search for second generation loptoquarks (LQ) produced in jfp collisions at 
^/s =1.8 TeV using the D0 detector at Fermilab. Second generation leptoquarks are assumed to 
be produced in pairs and to decay to either /j. or v and either a strange or a charm quark (q). Limits 
are placed on a(pp —> LQLQ ^ iiv + jets) as a function of the mass of the leptoquark. For equal 
branching ratios to ^q and vq, second generation scalar leptoquarks with a mass below 160 GeV/c^, 
vector leptoquarks with anomalous minimal vector couplings with a mass below 240 GcV/c^, and 
vector leptoquarks with Yang- Mills couplings with a mass below 290 GeV/c^, are excluded at the 
95% confidence level. 



2 



Leptoquarks (LQ) are hypothetical particles that carry 
color, fractional electric charge, and both lepton and 
baryon number. They appear in several extended gauge 
theories and composite models beyond the standard 
model 0]. Leptoquarks with universal couplings to all 
lepton flavors would give rise to flavor-changing neutral 
currents, and are therefore tightly constrained by exper- 
imental data To satisfy experimental constraints on 
flavor-changing neutral currents, leptoquarks that couple 
only to second generation leptons and quarks are consid- 
ered. 

This Letter reports on a search for second genera- 
tion leptoquark pairs produced in pp interactions at 
a center-of-mass energy y/s = 1.8 TcV. They are as- 
sumed Q to be produced dominantly via the strong 
interaction, pp ^ g + X LQLQ + X. The search is 
conducted for the signature where one of the lepto- 
quarks decays via LQ — )■ muon + quark and the other 
via LQ — > neutrino + quark, where the quark may be ei- 
ther a strange or a charm quark. The corresponding ex- 
perimental cross section is 2/3(1 — f3) x a{pp LQLQ) 
with (3 the unknown branching fraction to a charged lep- 
ton (e, /i, r) and a quark (jet) and (1 — /3) the branching 
fraction to a neutrino (v) and a jet. The search con- 
siders leptoquarks with scalar or vector couplings in the 
fj^iy + jets final state. Additional details on this anal- 
ysis may be found in reference ^ Previous studies by 
the D0 [^,|^ and CDF collaborations have considered 
the /i/i 4- jets final state for scalar couplings, resulting 
in limits of 140 GeV/c^ and 160 GeV/c^ respectively for 
P = 1/2. 

The D0 detector Q consists of three major compo- 
nents: an inner detector for tracking charged particles, a 
uranium-liquid argon calorimeter for measuring electro- 
magnetic and hadronic showers, and a muon spectrom- 
eter consisting of a magnetized iron toroid and three 
layers of drift tubes. Jets are measured with an en- 
ergy resolution of approximately fj{E) = 0.8/ \/E {E in 
GeV). Muons are measured with a momentum resolution 
cr(l/p) = 0.18(p- 2)/p2© 0.003 {p in GeV/c). 

Event samples are obtained from triggers requiring 
the presence of a muon candidate with transverse mo- 
mentum pt^ > 5 GeV/c in the fiducial region jry^l < 1.7 
{rj = — ln[tan(i0)], where 9 is the polar angle of the track 
with respect to the z axis taken along the proton beam 
line), and at least one jet candidate with transverse en- 
ergy Ei^ > 8 GeV and \rij\ < 2.5. The data used for this 
analysis correspond to an integrated luminosity of 94±5 
pb^i collected during the 1993-1995 and 1996 Tevatron 
collider runs at Fermilab. 

In the final event sample, muon candidates are required 
to have a reconstructed track originating from the inter- 
action region consistent with a muon of pj. > 25 GeV/c 
and \ri^\ < 0.95. To reduce backgrounds from heavy 
quark production, muons must be isolated from jets 
( A7^(^, jet) > 0.5 for E^. > 15 GeV, where A7^(^, jet) 
is the separation between the muon and jet in the rj — (p 



plane), and have energy deposition in the calorimeter 
consistent with that of a minimum ionizing particle. 
Events are required to have one muon satisfying these 
requirements. Events containing a second muon which 
satisfy these requirements, with the fiducial requirement 
relaxed to [yy^l < 1.7, are rejected. 

Jets are measured in the calorimeters and are recon- 
structed using a cone algorithm with a radius TZ = 0.5 
{JZ = _|_ Ary^). Jets must be produced within 

\rjj\ < 2.0, and have Ei^p > 15 GeV; with the most ener- 
getic jet in each event required to have |77j| < 1.5. 

The transverse energy of the neutrino is not directly 
measured, but is inferred from the energy imbalance in 
the calorimeters and the momentum of the reconstructed 
muon. Events are required to have missing transverse 
energy > 30 GeV. To ensure that frp is not dominated 
by mismeasurement of the muon px, events having ]/Jj^ 
within TT ± 0.1 radians of the muon track in azimuth are 
rejected. 

To provide further rejection against dimuon events in 
which one of the muons was not identified in the spec- 
trometer, muons are identified by a pattern of isolated 
energy deposited in the longitudinal segments of the 
hadronic calorimeter ^ . Any event where such deposited 
energy lies along a track originating from the interaction 
vertex in the region \r]\ < 1.7 and is within 0.25 radians 
in azimuth of the direction of the vector is rejected. 

Each candidate event is required to pass a selection 
based on the expected LQ event topology. Since the de- 
cay products of the LQ are /iq or vq, the muon and neu- 
trino in LQ pair decays come from different parent parti- 
cles nearly at rest and are therefore uncorrelated. For the 
primary background events (e.g. W + jets), the two lep- 
tons have the same parent. Similar reasoning holds for 
the jets. Correlated backgrounds are rejected with the 
requirement of significant separation between the muon 
and ]/]rp (|A0(/i, $rp)\ > 0.3) and between the two leading 
jets (A7^(Jl,j2)> 1.4). 

The ISAJET ||l^ Monte Carlo event generator is used to 
simulate the scalar leptoquark {Slq) signal, and PYTHIA 
jm is used for the vector leptoquark (Vlq) signal. The 
efficiencies for Vlq and Slq are consistent within dif- 
ferences due to the choice of generator. This is verified 
by choosing a test point at which both scalar and vector 
Monte Carlo events from the same generator are com- 
pared. Therefore, efficiencies obtained from the two sim- 
ulations are not distinguished. In addition, the efficien- 
cies for vector leptoquarks are insensitive to differences 
between minimal vector (kg 1; Aq = [^) and Yang- 
Mills {kg — 0; Xg — [|2|) couplings at large mass Q 
{Mvj^Q > 200 GeV/c^). The leptoquark production cross 
sections used for the Slq are from next-to-leading or- 
der (NLO) calculations jlSj with a renormalization scale 
fi = Ms^Q and uncertainties determined from variation 
of the renormalization/factorization scales from 2Ms^q 
to iMsj^Q. The Vlq cross sections are leading order (LQ) 
calculations at a scale /i — Mvj^q [Ol. 
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FIG. 1. Kinematic distributions for pjy + jets events. The 
quantities shown in (a)-(d) are used as inputs into the neural 
network (see text). The shaded regions give the background 
expectations, the square points are the + jets data, and 
the triangular points are signal Monte Carlo. 

The dominant backgrounds, from W + jets and 
Z + jets, are simulated using VECBOS for parton 
level generation and HERWIG for parton fragmen- 
tation. Background due to WW production is simulated 
with PYTHIA [|ll| . Additional background from tt decays 
into one or more muons and two or more jets, is simu- 
lated using the HERWIG Monte Carlo program for a top 
quark mass of 170 GeV/c^. Monte Carlo samples are 
processed through a detector simulation program based 
on the GEANT package. 

With the initial data selection described above, there 
are 107 events, consistent with a background of 106±30 
events (see Fig. 1). The dominant background is 
W + jets with 100±30 events. Other backgrounds are 
2.7±0.7 {Z + jets), 2.4±0.8 ti, and 1.5±0.6 (WW). The 
uncertainty in the background is dominated by the sta- 
tistical uncertainty in the W + jets simulation and the 
systematic uncertainty in the W + jets cross section. 
The expected signal for 160 GeV/c^ scalar leptoquarks is 
4.8±0.7 events. Signal estimations are shown for a Slq 
mass of 160 GeV/c^ using the NLO cross section with a 
scale of 2Ms^Q- 

To separate any possible signal from the backgrounds, 
a neural network (NN) |l^ with inputs: Eij^ , E'ij^ , and 
Iprp and nine nodes in a single hidden layer is used. The 
network is trained on a mixture of -I- jets, Z -\- jets 
and tt background Monte Carlo events, and an indepen- 
dently generated signal Monte Carlo sample at a mass 
of 160 GeV/c^. Figure 1 shows distributions of the four 
input quantities and Fig. 2 the network output (referred 
to as the discriminant, I?nn)- No evidence of a signal is 
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FIG. 2. Output of the neural network. The network calcu- 
lates a value for each event based on the inputs (shown in Fig. 
1) and a set of internal values which are determined during 
network training on signal and background Monte Carlo. 

observed in either the discriminant distribution or any of 
the kinematic distributions. For setting limits, the selec- 
tion on Dnn is optimized by maximizing a measure of 
sensitivity ||T^ defined by 

n 

S{Dn^) = P{k, b)Mf^'{k, b, s{Mlq)) 

fe=0 

where P{k,b) = e-^b^/k\ is a Poisson coefficient with k 
being any possible number of observable events, b the ex- 
pected mean number of background events, and s{Mlq) 
the expected signal for a given leptoquark mass. 

^95% 

is an approximate mass limit at the 95% confidence 
level for a given fc, s and b. S'(il'NN) is the sum of the 
approximate mass limits, weighted by the probability of 
observing fc = 0, 1, 2, . . . , n{P{n, b) < 0.05) events for a 
particular choice of the Dnn selection criterion. 

By maximizing the value of S'(Dnn) a- discriminant 
selection of Unn > 0.9 is obtained. With this selection, 
no events remain in the data, which is consistent with an 
expected background of 0.7±0.9 events. The remaining 
background is dominated by ti (0.6±0.2 events). The 
uncertainty on the total background is dominated by the 
statistical and systematic uncertainties from W + jets. 

Table I shows the signal detection efficiencies and up- 
per limits ll^] on the cross section at the 95% confidence 
level as a function of the leptoquark mass. The domi- 
nant systematic uncertainty on the signal efficiency is due 
to the simulation, (initial and final state radiation, par- 
ton distribution function, renormalization scale, choice 
of generator) with a 10% uncertainty. The systematic 
uncertainties shown include approximately equal contri- 
butions from uncertainty in the jet energy scale |2l| and 
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LQ Mass 


efficiency 


^95% 




BRXUmV 


BRXCTy 


(GeV/c^) 


(%) 


(Pb) 


(pb) 


(pb) 


(pb) 


100 


3.7±0.2±0.6 


0.94 


2.8 


53 


430 


120 


5.0±0.2±0.7 


0.72 


2.2 


23 


150 


140 


7.2±0.3±1.1 


0.47 


0.75 


10 


50 


160 


10.3±0.3±1.5 


0.33 


0.34 


4.0 


25 


180 


12.2±0.3±1.8 


0.27 


0.16 
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10 


200 


13.4±0.3±2.0 


0.25 


0.08 


1.0 
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220 


14.1±0.3±2.1 


0.24 


0.04 


0.45 
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240 


15.2±0.3±2.3 


0.23 


0.02 


0.23 


1.3 
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15.5±0.3±2.3 


0.22 


0.01 


0.13 


0.60 


280 


16.3±0.4±2.4 


0.21 




0.06 


0.30 


300 


15.7±0.4±2.3 


0.22 




0.03 


0.18 


350 


16.4±0.4±2.4 


0.21 






0.03 


400 


17.2±0.4±2.6 


0.20 









TABLE I. Signal detection efficiencies (with statistical and 
systematic uncertainty) and cross section limits (95% CL) for 
leptoquarks in the fiv + jets decay channel. Also shown for 
comparison are the expected cross sections for (3 — ^. (Js,;^q 
denotes the theoretical cross section for scalar leptoquarks 
with a scale 2Msi^q , (Tmv the cross section for vector lepto- 
quarks with anomalous minimal vector couplings, and ayM 
leptoquarks with vector Yang-Mills couplings. 

the trigger efRciency/spectrometer resolution for high pt 
muons (6.6% and 6.4% respectively). The overall sys- 
tematic uncertainty for the signal elRciency is 15%. 

The limits on the observed cross section are shown in 
Fig. 3, and are compared with the theoretical cross sec- 
tion times branching ratio for scalar and vector lepto- 
quark production for /? = i. Mass limits of 160 GeV/c^ 
for scalar leptoquarks and 290 (240) GeV/c^ for vector 
leptoquarks with Yang-Mills (minimal vector) couplings, 
are obtained at the 95% confidence level. 

In conclusion, we have performed a search for second 
generation leptoquarks in the ^u+jets decay channel us- 
ing 94 ±5 pb"-'^ of data collected with the D0 detector at 
the Fermilab Tevatron. No evidence for a signal is seen 
and limits are set at the 95% confidence level on the mass 
of second generation leptoquarks. For equal branching 
fractions to fj,q and jyq (/3 = i) hmits of 160 GeV/c^, 
240 GeV/c^, and 290 GeV/c^ for Slq, minimal vec- 
tor, and Yang-Mills vector couplings, respectively, are 
obtained. 
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FIG. 3. Cross section limits in the fiv + jets channel. The 
Vlq cross sections are leading order[12], calculated at a scale 
jj, = Mvj^Q ■ The Slq cross sections are next-to-leading or- 
der [13]. The calculation is done at a ronormalization scale 
jj, = Mslq with uncertainties obtained from variation of the 
renormalization/factorization scale from 2Msj^q to ^Ms^^. 
For the Slq the limit is obtained at the intersection of the 
experimental curve with the theoretical curve for fx = 2Ms^g . 
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